Abnormal Electronic Transport in Disordered Graphene Nanoribbon 
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We investigate tlie conductivity a of graphene nanoribbons with zigzag edges as a function of Fermi 
energy Ef in the presence of the impurities with different potential range. The dependence of o{Ef) 
displays four different types of behavior, classified to difi'erent regimes of length scales decided by the 
impurity potential range and its density. Particularly, low density of long range impurities results in 
an extremely low conductance compared to the ballistic value, a linear dependence of a{Ep) and a 
wide dip near the Dirac point, due to the special properties of long range potential and edge states. 
These behaviors agree well with the results from a recent experiment by Miao et al. (to appear in 
Science) . 
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Introduction. — Recent breakthrough in graphene fab- 
rication has attracted many attentions to this two- 
dimensional (2D) material [l]. The honeycomb lattice 
structure of graphene gives rise to two interesting elec- 
tronic properties in the low energy region which are 
distinct from conventional 2D materials, i.e., two val- 
leys associated to two inequivalent points K and K' 
at the corner of the Brillouin zone, and linear "Dirac- 
like" rather than quadratic bare kinetic energy dispersion 
spectra. Many of the interesting experimental results are 
attributed to these peculiar properties near Dirac point, 
the Fermi level for undoped graphene. 



Some interesting aspects of the electronic transport in 
disordered graphene have been investigated theoretically 

0, [3, 0, [1| and experimentally [1, 0, Bl^l • It was realized 
that the potential range of the impurities plays a spe - 
cial role in the electronic transport in graphene 1^, 11|. 
Impurities with long range potential scattering were con- 
sidered to be a possible origin for some unconventional 
features in the experiments ilS, J3, 14, 15, 16, J, 7]. Such 
a potential could be realized by screened charges in the 
substrate. The peculiarity of the long-range disorder is 
the absence of valley mixing due to the lack of scatter- 
ing with large momentum transfer. In a realistic experi- 
ment, a gate voltage Vg can continuously tune the carrier 
density (thus the Fermi energy Ep) in the graphene sam- 
ple. A perfect linear relation between conductivity a and 
gate voltage Vg was observed [l[. However, clear nonlin- 
ear cr{Vg) curves emerge in a recent experiment 18|. For 
large Vg, the cr(V^) curves show a sub-linear behavior, 

1. e., square root in Vg rather than a linear one. The con- 
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ductance is smaller than the theoretical ballistic value by 
a factor of 3-10. Whereas in the low Vg region near the 
Dirac point, even this square-root like behavior breaks 
down, and a wide dip appears. This dip is wider for a 
smaller sample. These observed novel transport features 
have no explanations thus far. 




FIG. 1: (color online) Geometry of the device for measuring 
the conductance of graphene nanoribbon. Two semi-infinite 
leads (red) are connected to the graphene nanoribbon sam- 
ple with zigzag edges (black). The coordinate of each site is 
labelled as {ix,iy), where 1 < < and 1 < iy < Ny are 
integers. 



In this Letter, we perform systematic calculations to 
investigate the effect of the impurity potential range and 
its density on the conductance of graphene nanoribbon. 
The dependence of a{Ep) displays four different types of 
behavior, corresponding to regimes with different length 
scales depending on the range and density of the impuri- 
ties, which can be used as criterions in the experiments. 
Moreover, we demonstrate that the nonlinear dependence 
in the recent experiment [l8l| can also be explained when 
scattering due to the low density and long range impuri- 
ties are accounted. 

Model and Method. — We consider a two-terminal de- 
vice to calculate the conductance, which includes a 
graphene nanoribbon and two leads with zigzag edges, 
as shown in Fig. [U where the graphene sample is di- 
vided into vertical chains with Ny sites in each chains. 
In this setting, the coordinate (ix, iy) (1 < ix < N^, 
1 < iy < Ny) of each site is labelled. Two clean and 



FIG. 2: (color online) (a) Conductivity a of graphene nanorib- 
bons with zigzag edges as a function of Fermi energy Ef with 
different potential range 5. (b) Conductivity cr as a function 
of potential rnge ^ at different Fermi energy Ef- Nx ~ 106, 
Ny = 60 {Lx = 90.9a and Ly = 89.0a), W = 0.5t for both 
cases, where t = 2.7eV and a = 1.42nm. The conductivity is 
averaged over 100 random configurations for each curve. 



semi-infinite leads are assumed to have the same type of 
lattice as the graphene sample 1^, 2^ to avoid additional 
scattering contribution from the mismatched interfaces 
between different types of lattices. 

We describe graphene by the tight binding Hamilto- 
nian for the tt orbital of carbon 12 ll. 12211 



i ii,3) 



(1) 



where cj (c^) creates (annihilates) an electron on site i, 
Si is the potential energy and t (~2.7eV) is the hopping 
integral between the nearest neighbor carbon atoms with 
distance a ('^1.42 A). We use t as the energy unit and a 
as the length unit. 

In the presence of disorder, Ni impurities are randomly 
distributed among N (= x Ny) sites. The potential 
energy of the i-th site £i at position is induced by these 
impurities as [l^. \y\ 



Ni 

^Kexp(-|ri-r„|V(2a), 



(2) 



where rn is the position of the n-th impurity, ^ repre- 
sents the spatial range of the impurity potential, and the 
potential strength Vn of the impurities is randomly dis- 
tributed in the range (—14^/2, 1^/2), independently. The 
average distance between two impurities Ri can be de- 
fined as Ri ^ ^ LxLy/Ni, where ^^(y) is the length 
(width) of the rectangular sample. We shall show in the 
following that distinct interesting phenomena can be ob- 
served when the system is in the different regimes of these 
length scales that is determined by Lx(y), C ^^id Ri. 

In the framework of the non-equilibrium Green's 
function method , the zero temperature conduc- 
tance G and density of states (DOS) p of the 



FIG. 3: Conductivity of graphene nanoribbons with zigzag 
edges as a function of Fermi energy, with A'^^, = 106, Ny — 60 
{Lx = 90.9a and Ly — 89a), W — 0.5t at different ^ and Ui. 
The conductivity is averaged over 100 random configurations 
for each curve. Note different a scalings in each regime, (a) 
Regime 1, no impurity, (b) Regime 2, short range impurities, 
(c) Regime 3, long range and low density impurities, (d) 
Regime 4, long range and high density impurities. 



sample at Fermi energy Ep can be written as 
GiEp) = ^TTilm^l{EF)G-{EF)linJ:'-j,{EF)G%EF)) 
and p{Ef) = -ilm[TrG"^(£;F)], where G''^''\Ef) is the 
retarded (advanced) Green's function, and S^^^"^) is the 
retarded (advanced) self-energy due to the left (right) 
lead O, [ii, [23, [ii]. The conductivity is related to the 
conductance by the geometric relation a — L^GjLy. 

In the clean limit, a self-consistent calculation for 
graphene in the Hartree approximation shows that, a oc 
^y [ii]. In this case, a Ef (also see Fig. [3] 



(a)). This leads Xq Ef — ay 



where ay is a device- 



dependent prefactor, whose typical value ~ 10^ — 10^ . 
This relation between Ef and Vg is valid even in the pres- 
ence of disorder since it is a global response. Therefore, 
we can concentrate on the relation between a and Ef. 

Results and Discussions. — Firstly, we investigate the 
effect of the potential range ^ and density Ui = Ni/N 
of the impurities with fixed W = 0.5t. Let us start by 
having a first glance at the effect of ^. In Fig. [2] (a), 
we plot a{EF) with different ^. A direct conclusion from 
this figure is that, in most regions of Ef, short range 
(^ = 0, la) impurities can lead to a considerable decrease 
of conductivity, but long range (^ = 7a, 15a) impurities 
will decrease the conductivity much further, in the case 
of same impurity density. In the experiments by Miao et 
al. the conductance is smaller than its ballistic value by a 
factor of 3-10, in the range of high Vg [ll]. This suggests 
that the samples used in these experiments must include 
long range impurities. Near the Dirac point Ef = 0, 
this rapid variation breaks down, see Fig. [2] (b). The 
magnitude of <j{Ef = 0) is relatively universal /h), 
compared to a at finite Ef. This is consistent with the 
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FIG. 4: (color online) (a) Averaged conductivity of graphene 
ribbons with short range impurities (black) (Data are same 
with Fig[3](b)) and long range impurities (red) (Data are same 
with Fig|3] (c)), and density of states of a clean ribbon with 
the same size (green, in arbitrary unit) . (b) Band structure of 
a zigzag edge graphene nanoribbon with the same size. Note 
correspondences between these figures. 
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FIG. 5: (color online) Conductivity ct as a function of Vg: A 
qualitative comparison between theory and experiment. The 
black lines represent experimental results by Miao et al. p^ . 
The upper curve is for a larger sample and the lower curve is 
for a smaller sample. The red lines represent our numerical 
results. The upper curve: A^^: = 120, Ny = 110, ^ = 4a, 
W = 2t; The lower curve: A^^ = 72, Ny = 80, C = 4a, W = 3t. 
The conductivity is averaged over 50 random configurations 
for each curve. 



experiment fl]. The linear dep endence of u{Ef) obtained 
in the mean field theory [28| is not entirely vahd. Our 
numerical results indicate that a{Ep) is a more compli- 
cated function that depends on the nature of disorder. 
To illustrate the physics, we focus on the conductivity a 
for graphene nanoribbons with zigzag edges for different 
^ and rii, which can be seen in Fig. [3l All these behaviors 
can be classified into four typical regimes. 

Regime 1, no impurity, rii = 0. When there is no impu- 
rities {rii — 0), <7{Ep) increases almost linearly except for 
small quantized plateaus due to finite size quantization 
in the transverse direction 27|, as shown in Fig. [3] (a). 
When the sample is large enough, these sub-structures 
due to finite size effect can be ignored, so ajErr) is lin- 
ear, and a ~ ^/Vg according to Ep — \/^ [26|. 

Regime 2, short range impurities, ^ a. When short 
range impurities are present (Fig. [3](b)), the first distinct 
feature is the sub-linear behavior of a(Ep), especially 
in the high energy region. This anomaly can be under- 
stood as enhanced scattering due to large level broaden- 
ing when DOS is large. In the coherent phase approxi- 
mation without vertex corrections, the conductivity a of 
disordered system can be written as [1^ 



bp e vp 



(3) 



where Sp is the area of the Fermi sphere and S2 is the 
imaginary part of self energy corresponding to the disor- 
der scattering. The level broadening IS2I = TrpNi\V{q)\'^ , 
where V{q) is the Fourier transform of impurity potential 
[2^. For a disordered nanoribbon, the van Hove singu- 
larity at E = t oi 2D graphene [2^, [s^l degenerates into 



a finite but still very sharp peak (see the green curve in 
FigHl (a)). Therefore, E2 also has a sharp peak at this 
point, giving rise to a minimum of a as can be seen from 
This minimum a at Ep = t leads to a sub- linear 
a{Ep) and sub-square root o^{Vg), according to the rela- 
tion between Ep and Vg mentioned above. 

Such level-broadening enhanced scattering also hap- 
pens at the bottoms of sub-bands, where van Hove sin- 
gularities emerge [3^ 31 1. Indeed, when disorder is 



not strong enough to smear these singularities out com- 
pletely, a small dip can be observed at each sub-band 
bottom, as can be seen from Fig. |4] (a). 

Variation of of short range impurities does not 
change the qualitative behavior of a{Ep), but reduces 
the magnitude of a for a given Ep, when ^ 1, as one 
expects. However, for extremely short range impurities 
{£, ~ 0), we find that even the magnitude of (j{Ep) is 
quite independent of Ui when Ui > 20%. 

Regime 3. long range and low density impurities, ^ 3> 
a and ^ ^ Ri. When the potential range ^ increases 
further, interesting physics appears. As can be seen in 
Fig. [3](c), the a{Ep) curves resume their linear behavior 
in most energy regions (while the slope is much smaller 
as mentioned above). This manifests suppression of large 
momentum scattering due to the long range impurities. 

Another notable nonlinear a{Ep) can be observed near 
the Dirac point, where a wide dip appears. We find 
this happens within the energy region where the first 
sub-band is visible for a clean graphene (see Fig[4] (b)). 
For a smaller graphene sample, different sub-bands are 
more separated in energy than for a larger sample, giving 
rise to wider quantized conductance plateaus, and also a 
wider dip near the Dirac point, which is consistent with 
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the experimental result [l8|. 

This wide dip can be understood as follows. As shown 
in Fig. m (b), the band structure near the Dirac point 
is composed of two branches of subbands, i.e., the up- 
per band +E{kx) > and the lower band —E{kx) < 0. 
These branches correspond to binding and antibinding 
states localized at different edges and sublattices [H, . 
The wave functions of these edge states possess different 
signs according to two edges, sublattices and branches 
34| . Long range impurities will (while short range ones 
will not) couple these edge states, giving rise to rather 
large scattering matrix elements between two valleys. 
When E = 0, two branches and valleys degenerate, the 
magnitude of scattering matrix elements decreases since 
different signs of these edge states at E = making a 
larger possibility of canceling each other. This is verified 
by numerical calculations for scattering matrix elements. 

Therefore, ^{Vg) behaves in a square root way in this 
regime, by noting Ep — ctVg \f^g- In Fig- El we plot 
ct(V^), setting avg = 10~^. A perfect quantitative fitting 
cannot be reached because the size of the sample in the 
experiments is the order of Nx.y ~ 10'^, which is well 
out of the capability of numerical calculations. But the 
qualitative features, i.e., sub-linearity and wide dip in the 
experiment (Fig. [5]) can be clearly seen. Once again, we 
attribute this experimental result to the contribution of 
the low-density and long-range impurities. 

Regime 4, long range and high density impurities, ^ ^ 
a and ^ ^ Ri. In the case of low density impurities, the 
scatterings due to different impurities are independent. 
But when the density is sufficiently high, so that po- 
tential field induced by different impurities overlap, and 
multi-scattering dominates. This multi-scatterings have 
no obvious effect on the existence of the dip. While in 
the high energy region, the linear relation <t{Ep) breaks 
down and the curve degenerates into a square root like 
curve and cr(Vg) ~ ^/Vg correspondingly, see Fig. [3] (d). 

Finally, we discuss the energy scaling of the impurity 
potential W considered to be fixed thus far. In our cal- 
culations W is much larger than the level spacing of sub- 
bands. The opposite limit has been investigated recently, 
and a perfectly conducting channel was found [l7| . 

Conclusions. — As a summary, we numerically investi- 
gate the transport properties of graphene nanoribbons 
in the presence of the impurities with different density 
and potential range. In the Fermi energy region of fo- 
cus, four typical types of behavior can appear from the 
unconventional electronic structures in zigzag graphene 
nanoribbons, which can be tested by future experiments. 
The third regime for the low density and long range im- 
purities can be used to exp lain the nonlinearity of <j{Vg) 
in a recent experiment [18|. 
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